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Abstract 



The usual dispersionless limit of the KP hierarchy does not work in the case where the dependent 
' variable has values in a noncommutative (e.g. matrix) algebra. Passing over to the potential KP 

. hierarchy, there is a corresponding scaling limit in the noncommutative case, which turns out to 

' be the hierarchy of a "pseudodual chiral model" in 2 + 1 dimensions ("pseudodual" to a hierarchy 

. extending Ward's (modified) integrable chiral model). Applying the scaling procedure to a method 

generating exact solutions of a matrix (potential) KP hierarchy from solutions of a matrix linear heat 
\^ . hierarchy, leads to a corresponding method that generates exact solutions of the matrix dispersionless 

' potential KP hierarchy, i.e. the pseudodual chiral model hierarchy. We use this result to construct 

. classes of exact solutions of the su{m) pseudodual chiral model in 2 + 1 dimensions, including 

various multiple lump configurations. 



1 Introduction 

Expressing the scalar KP hierarchy with dependent variable u{ti ,t2, ■ ■ terms of new evolution vari- 
ables Tn = etn with a parameter e, the limit e ^ (keeping T„ fixed) leads to the so-called disper- 
sionless KP hierarchy (see [1-23], for example). The same hmit does not work, however, for the KP 
hierarchy with dependent variable in a noncommutative (e.g. matrix) algebra. In fact, different scaling 
limits of the matrix KP equation have already been explored in [24], where the multiscale expansion 
method has been used to relate different integrable systems. 

In the present work we formulate a dispersionless limit of the "noncommutative" potential KP 
(ncpKP) hierarchy with dependent variable 0, where u = 4>ti- It turns out to be the hierarchy asso- 
ciated with a "pseudodual chiral model" (pdCM) in 2 + 1 dimensions, a well-known reduction of the 
self-dual Yang-Mills equation [25, 26]. Applying the scaling limit procedure to a method generating 
exact solutions of a matrix pKP hierarchy from solutions of a matrix linear heat hierarchy, then results in 
a method generating solutions of this pdCM hierarchy. 
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In section [2] we consider the dispersionless limit of the ncpKP equation. Section [3] generalizes 
this limit to the whole ncpKP hierarchy, explores some of its properties, and in particular establishes 
a pseudoduality relation with a hierarchy that extends Ward's (modified) chiral model in 2 + 1 dimen- 
sions [27-32]. The latter model has been studied extensively [33-52] (see also [53-55] for the Ward 
model in (anti-) de Sitter space-time and [56-63] for explorations of a Moyal-deformed version), in 
particular concerning its (multi-) lump solutions, which are two-dimensional soliton-like objects. In 
this respect, its pseudodual received comparatively little attention. The dependent variables of the two 
equations are related by a kind of hetero-Backlund transformation. Given a solution of one of the two 
equation, this becomes a first order system of partial differential equations, which determines a solution 
of the other equation. The necessary integration is typically difficult to carry out, however. Hence, al- 
though some properties of the pseudodual model can certainly be infered from corresponding knowledge 
of the Ward model, there is no explicit translation of its solutions. In any case, in this work we present 
an independent approach to solutions of the pdCM and moreover to its hierarchy. 

In section |4] we derive the abovementioned method to generate exact solutions of the pdCM hierar- 
chy from corresponding knowledge of the ncpKP hierarchy. The main result is independently verified in 
section |5] and then applied to construct some classes of exact solutions. This section is actually formu- 
lated in such a way that it can be accessed almost without any knowledge of the previous sections. We 
concentrate on solutions of the su{m) pdCM hierarchy and restrict concrete examples to the su{2) case. 
Some conclusions are collected in section |6] 

2 The dispersionless limit of the noncommutative pKP equation 

Let (f){t) with t = {ti,t2, . . .) be a function with values in some matrix spaceQ .4 which is endowed with 
a product A - B = AQB, where Q is a constant matrix, i.e. independent of t. We consider the following 
ncpKP equation, 

^(t>tx - (f>xxxx -^(t>yy = Q{4>xQ(t>x)x - 6 [4>x,(t>y]Q , (2-1) 

where x = ti,y = t2,t = and 

[A,B]q:=AQB-BQA. (2.2) 
Let (p now also depend on a parameter e in such a way that 

0(t,e) = e"«>(T) + 0(e"+^) (2.3) 
with some integer a. Furthermore, we assume that Q has an expansion 

Q = Q(0) + eQ(i) + O(e') • (2.4) 

Rewriting the ncpKP equation in terms of the rescaled variables Tn = etn, dividing the equation by the 
maximal power of e common to all of its summands, and taking the limit e ^ while keeping Ti,T2, . . . 
fixed, should result in an equation that still has linear as well as nonlinear terms (in $). This fixes the 
value of a, but we have to distinguish the following two cases. 

If the algebra {A, •) is commutative at e = with (5(o) / 0, and hence the commutator , ^y]q^q^ = 
^xQ{o)^Y — ^yQ{o)^x vanishes, then our requirements lead to a = — 1, and the scaling limit of the 
pKP equation, divided by e, is 

A^TX - 3$yy = 6($xQ(o)^x)x - 6 [<^x ,'^y]q^,^ , (2.5) 

'The entries will be taken as complex functions of ti, i2, • ■ ■, though large parts of this work also apply to the case where 
they are elements of any (possibly noncommutative) associative algebra, for which differentiability with respect to ti,t2, ■ ■ ■ 
can be defined. 
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where X = Ti,Y = T2,T = T3. If ^ is a scalar and (5(o) = 1, the last equation reduces to 

4$TX-3«>yy = 6($x^)x . (2.6) 

This is the potential form of the dispersionless limit of the ("commutative") scalar KP equation, which 
is also known as the Khokhlov-Zabolotskaya equation (see [4] for instance). 
If the algebra {A, •) is noncommutative at e = 0, we have to set 

a = , (2.7) 

and this choice will be made throughout this work. Then we obtain the following dispersionless limit of 
the ncpKP equation (12. IK 

4$TX -3$yy = -el^X.-^ylQjo) ■ (2.8) 

Up to the modified matrix product and rescalings of the coordinates, this is a well-known reduction of the 
self-dual Yang-Mills equation (see [25,26,64-68]). With the further dimensional reduction <^x = ^T> 
it becomes the pseudodual chiral model [69-71] (see also [64,72-74]). Accordingly, we may call (12.81) a 
pseudodual chiral model in 2 + 1 dimensions, in the following abbreviated to pdCM. In fact, as explained 
in section I3l2l it is "pseudodual" to an integrable (modified) chiral model in 2 + 1 dimensions. 



3 The dispersionless limit of the ncpKP hierarchy 

A functional representation of the ncpKP hierarchy is given by [75] 

- 0-[A])(A-^ - Q4>) -<\>t,=B- B-xx] , (3.9) 

where Q is an arbitrary ,A-valued function, and (0_[A])(t) := (/)(t — [A]) is a Miwa shift with [A] = 
(A, A^/2, A'^/S, . . .), A an indeterminate. Eliminating Q from this equation, we get the following func- 
tional form of the ncpKP hierarchy, 

((0 - 0-[A])(A-^ - Q</.) -4>t^- ((0 - 0-[A])(A-^ - Q^) - 

= ((0 -0-m)(m"'-Q</') ((<^-</'-M)(/^"'-W)-'^ti)_j^j> (3-10) 

where /i is another indeterminate. 

If p„, n = 1, 2, . . ., denote the elementary Schur polynomials, then 

p,(-a) = --aT„+0(e2), (3.11) 
n 

where d = {dt^ , dt^l2, dt.^/2>, . . .), and hence 

0-0_[A] = eZ?(A)cI> + 0(e2), (3.12) 

where 

^(^)^=E-^^- (3.13) 

n>l 

In accordance with (12.31) . where now a = 0, we shall assume 

0(t,e) = G(T) + 0(e) . (3.14) 
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Then we obtain 

D{X){<!>) (A-i - Q(o)«>) - «>ri = ^(A)(e) . (3.15) 
Expanding this in powers of A, we find 

^^T.^i - ^^T„Q(o)$ = ^eT„ n = 1, 2 . . . . (3.16) 
Elimination of © results in the hierarchy equations 

Ti TJX 

^^^:^^r„+i,T™ - (^ + i) ^W,T„ = ^T„Q(0)^T™ - ^T^Q(0)^T„ . (3.17) 

Introducing 

Xn--=nTn n = l,2, (3.18) 

this becomes 

= [^x„,^x„](3(o) m,n = l,2,.... (3.19) 

For m = 1, n = 2, we recover (12.81) . 
Expressing (5(o) 

Q(o) = yC/t (3.20) 
with matrices U, V and the adjoint (complex conjugate and transpose) C/^ of U, then 

:= C/t$y (3.21) 

(which includes the cases Lp = Q(q^^ and Lp = $(5(0)) solves 

'fxr. + UXm- 'Pxm + l,Xn = VPxr,,^xJ\ 771,71=1,2,..., (3.22) 

if $ solves (13.191 ). The power of this observation lies in the fact that any solution of (13.191 ) in some 
M X N matrix algebra, where (5(o) = VU^ with an M x m matrix U and an x m matrix V, 
determines in this way a solution of (13.221 ) in the m x m matrix algebra. For example, if we are looking 
for solutions of (13.221 ) in the algebra of 2 x 2 matrices, we may first look for solutions of (13.191 ) with any 
M,N >2 and (5(o) = VU"^ with M x 2 and N x 2 matrices U and V. In this way (simple) solutions of 
(13.191 ) in arbitrarily large matrix algebras lead to (complicated) solutions of (13.221 ) in the algebra of 2 x 2 
matrices. In particular, this explains the significance of Q^q) our previous formulae. In section |5] we 
will substantiate this method. The hierarchy (13.221 ) is consistent with restricting ip to take values in any 
Lie algebra, e.g. sl{N, M), sl{N, C), u{N) or su{N). If ip solves (13. 221 ). then also ip + ipQ, where tpo is 
a constant in the respective Lie algebra. 

As a consequence of their origin, the hierarchies ( 13.191 ) and (13.221 ) are invariant under the scaling 
transformation x„ A x„, n = 1, 2, . . ., with any constant A / 0. 

Remark. If 51, 52 are any two constant invertible matrices with size such that gi^g2 is defined, then 

^ ^ 91^92 , Q(o) ^ 92^Q{0)9i^ (3-23) 
leaves (13.191 ) invariant. If (5(o) is given by (I3.20I ). the latter transformation results from 

V^g^^V, U^iglr^U, (3.24) 

and ip is invariant. More generally, the transformation V g2^Va, U ^ {g\)^^U{a^)^^, with a 
constant m x m matrix a, leads to ip ^ a~^(pa. This leaves the hierarchy equations (13.221 ) invariant. □ 
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3.1 Some properties of the first dispersionless hierarchy equation 

A Lagrangian for the first equation (m = f , n = 2) 

^Xl,X3 - V>X2,X2 = - [^Xi,^X2] (3.25) 

of the hierarchy (13.221 ) is 

C = -tr(^v9a;iV7^3 - (/3^2^ - '^ip[(p^^,ip^2]) (3-26) 
(see also [64, 65]). After passage to the new coordinates x, y, t given by 

xi = ^{t-x), X2 = y, X3 = ^{t + x), (3.27) 
equation (13.251 ) becomes 

V^tt - fxx - fyy + [ft - fxjfy] = , (3.28) 

and the Lagrangian takes the form 

where we introduced the components t]'^'^ (with respect to the coordinates (x'^) = {t,x,y)) of the 
Minkowski metric in 2 + 1 dimensions, the totally antisymmetric Levi-Civita pseudo-tensor with e^^'^ = 
1, and a constant covector Vp with components (1,1,0). As a consequence of the translational invariance 
of the Lagrangian, the energy-momentum tensor 



--tr{ird^fd,ip + - fVpePf^^d^ipd^if) , (3.29) 



provides us with the conserved densities 

1 / 2 \ 

= -tr(^(ft<fx - ^fifx, fy]^ , T^2 = -tT{ipt(fy) . (3.31) 

Then also 

£ = T%- T^i = -hi[{^t - fx? + fy^] (3.32) 

is a conserved density. For any non-zero anti-Hermitian matrix, the trace of the square of the matrix is 
real and negative. Hence £ provides us with a non-negative "energy" density in the case where f takes 
values in the Lie algebra u{m) of the unitary group. 

For any infinitesimal symmetry 5(p = (with a parameter q) of the Lagrangian, there is a 

conserved current 

J'--=t4J^i^), (3.33) 



^d{df,ip) da 

i.e., f?^J^ = 0. A symmetry of the above Lagrangian is given by 5f = [C, (/9]a with any constant 
(anti-Hermitian) matrix C. Hence 

= -tr {^(\ip,ft\ + \{f^fy - 'iWyf + fyf^)) (3-34) 

is a conserved density. 
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3.2 Relation with Ward's chiral model in 2 + 1 dimensions 

The hierarchy (13.221) is related to the hierarchy of an integrable (modified) chiral model in 2 + f dimen- 
sions. First we note that (13.221) is the integrability condition of the Unear system 

Jx„+i = -Jy^x„ n = l,2,... (3.35) 

with some invertible J. Rewriting this as 

= -J~^Jx„+i n = l,2,..., (3.36) 

we find that (13.221 ) is automatically satisfied and the integrability conditions now take the form 

{J-'Jx„^,)x^-{J-'Jx^^,)x„=0 m,n = l,2,... . (3.37) 

In conclusion, solutions J of (13.371) are in correspondence with solutions ip of (13.221 ) via (13.361) . This 
correspondence is of a nonlocal nature. In particular, given a solution ip of (13.221 ). (13.361 ) does not directly 
determine J^^J^j^. We first have to solve (13.351 ) for J in order to be able to calculate this expression. 

(13.371 ) is immediately recognized as the dispersionless limit of the noncommutative modified KP 
hierarchy (see equation (4.12) in [76]). 

For m = 1 and n = 2, (13.371) reads 

{J'^Jx,)x^-{J''Jx,)x,=0. (3.38) 

This equation apparently first appeared in [77,78]. It is a reduction of the self-dual Yang-Mills equation 
(see [26,65,77], for example). In terms of the coordinates x, y, t given by (13.271 ). it takes the form 

iJ-^Jt)t - {J-^Jx)x - {J'^Jy)y + [J'^Jx, J-^Jt] = , (3.39) 

or in tensor notation (using the summation convention) 

+ = 0, (3.40) 

where /i, = 0,1,2, {rj^^'^) = diag(l, — 1, — 1), and e^^'^ is antisymmetric with e'^^ = — e^'^ = 1 and 
zero otherwise. We note that the bivector e^'^ breaks Lorentz invariance in 2 + 1 dimensions. Using the 
Lorentz invariant Levi-Civita pseudo- tensor and the constant unit co vector Va with components (0, 0, 1), 
it can be expressed as e^'^ = Va^"^'^- Another integrable equation is obtained if we choose Va to be 
timelike [27,38,78]. (13.401 ) is Ward's (2 + 1) -dimensional generalization of the chiral (or sigma) model 
[27-32], see also [33-48,50-52,68]. J can be consistently restricted to any Lie group, e.g. SL{N,R), 
SL{N,C), U{N) or SU{N). 

Remark. According to (13.361 ) we have J^^ Jy = (fx — and J^^Jt + J^^Jx = —(fy, in terms of the 
variables x,y,t given by (13.271 ). Hence 

£ = fward - tr {J-^Jt J-^Jx) , (3.41) 

where 

fWard = -^tr {{J-^Jtf + {J-^Jxf + {J-^Jy?) (3.42) 

is the energy density of Ward's chiral model. The difference between <?ward and £ is not a local expres- 
sion in terms of 99. The appendix attempts to further clarify the relation between Ward's chiral model 
and the pdCM hierarchy (and yet another version of it). □ 
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3.3 An associated bidifferential calculus 

On the algebra Aof m x m matrices with entries depending smoothly on xi, X2, • • •, we introduce two 
linear maps d,db>|l 

dV' = X] ' = ^ i^x„+i dxn ■ (3.43) 

n>l n>l 

By use of the graded Leibniz rule they extend to a (bi-) differential graded algebra and satisfy 

d2 = = dd + dd = , (3.44) 
and hence we have a bidifferential calculus. Dressing d by setting 

= dil) - All) , (3.45) 

with a 1-form A = ^^^^^ A„_|_idx„, we find that d,D yields again a bidifferential calculus (D^ = 
dD + Dd = 0), iff 

d^ = , d^ = ^ A yl (3.46) 

(see also [81]). These equations cover Ward's chiral model hierarchy as well as its pseudodual, which is 
the dispersionless ncpKP hierarchy. Indeed, solving the first equation by setting 

A = d^, (3.47) 

the second reproduces the pdCM hierarchy 

ddv? = (^Lp^dLp . (3.48) 
Alternatively, solving the second of equations (13.461 ) by setting 

A = -J-^dJ, (3.49) 

we recover the hierarchy 

d(J-^dJ) = (3.50) 

associated with Ward's chiral model. The relation between both hierarchies is given by 

J-^dJ = -d^j (3.51) 

(which is (13.361 )). This may be regarded as a "Miura transformation". The linear system associated with 
the bidifferential calculus is 

D^-Ad^ = 0, (3.52) 
with a parameter A. Taking components of the differential forms, this reads 

- An+i - Xd^M = n = 1, 2, . . . . (3.53) 



^We note that A = TZ o A where TZ is the linear left ^-module map determined by TZ{dx„) = dxn-i for n > 1, and 
TZ{Ax-i) — 0. This makes contact with Frolicher-Nijenhuis theory [79], see also [80]. 
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The integrability conditions now have the form 

[5x„+i - An+i - Xd^„,d^^^^ - Am+i - Xd:rJ = . (3.54) 

Its multicomponent version (and with m,n ^2) appeared in [82] (see (2. 1), (2.2), and also the references 
therein). 

Nonlocal conserved currents are obtained in the following way [81]. Let d^o = 0. As a consequence 
of the bidifferential calculus structure, there we. Xn, n = 1,2, . . ., such that 

jn+i ■■= Dxn = -dXn+1 n = 0, 1, . . . (3.55) 

iteratively determines Xn, n = 1,2,.... For example, starting with xq = ^ (the unit matrix), we get 
ji = DX = —dip (using (13.471 )). hence xi = + a with da = 0, and thus also da = 0. In the second 
step we have j2 = 0(99 + a) = d(/j — d(/? (99 + a), and the construction of the next current requires the 
integration of dx2 = d^? (93 + a) — dip. The constant a actually turns out to be redundant and should be 
set to zero. 

A Backlund transformation is obtained from 

{d-\-^I)){I + X-^B) = {l + \-^B){d-\-^I)'), D' = d-^', (3.56) 

with an operator B (see [83]). Expanding in powers of A^^, we find 

[d, i3] = D - D' , t)B = BV)' . (3.57) 

Assuming B{ip) = Bif) with a matrix B, this means 

d{B) = A' - A, d{B) = AB- BA' . (3.58) 

Using (13.471 ) and solving the first of these equations by setting B = ip' — ip — a with da = 0, we obtain 
from the second 

d{p' — p) = dp {p' — p — a) — {p' — p — a) dp' , (3.59) 

a Backlund transformation of the pdCM hierarchy. Alternatively, using (13.491 ) and solving the second of 
equations (13.581 ) by setting B = —J^^JCJ' with d/C = 0, the first becomes 

J'-^dJ' - J-^dJ = d{J-^JCJ') , (3.60) 

a Backlund transformation of the (modified) chiral model hierarchy (see also [84] for the case of the 
chiral model on a two-dimensional space-time). 

If Bij leads from an ith to a jth solution, a permutability relation is given by 

B12 + B2A = Bi3 + B34 , B12 B24: = Bi3 B34 (3.61) 

(see [83]). This determines algebraically a forth solution from a given (first) solution and two Backlund 
descendants of it (with different parameters). 

4 Toward exact solutions of the dispersionless ncpKP hierarchy 

In this section we start with a result that determines a large class of exact solutions of an ncpKP hierarchy 
and use the scaling limit toward the dispersionless hierarchy in order to obtain from it a corresponding 
result that determines exact solutions of the latter, which is a pdCM hierarchy. Let us recall theorem 4. 1 
from [75]. 
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Theorem 1 Let {A, •) be the algebra of M x N matrices of functions oft with the product 

A-B = AQB , (4.62) 

where the ordinary matrix product is used on the right hand side, and Q is a constant N x M matrix. 
Let X be an invertible N x N matrix and y £ A, such that X, y solve the linear heat hierarchy (i.e. 
dt„ (X) = dl[ (X), n = 2, 3, . . ., and correspondingly for y) and satisfy 

Xt,=RX + Qy, (4.63) 

with a constant N x N matrix R. The pKP hierarchy in {A^ •) is then solved by 

^■=yX-^. (4.64) 

□ 

A functional representation of the heat hierarchy condition is 

X-\X - X,^^-^) = Xt, , (4.65) 

and correspondingly for 3^ (with an indeterminate A). The theorem provides us with a method to construct 
exact solutions of the ncpKP hierarchy in (A, •). The idea is now to take the dispersionless limit of (14.631 ) 
and (14.651 ). This should then result in conditions that determine exact solutions of the pdCM hierarchy 
in {A, •). However, assuming for X,y power series expansions in e with nonvanishing terms of zeroth 
order, this results in too restrictive conditions. The way out is to note that a "gauge transformation" 

X = XG, y = yG, (4.66) 

with an X matrix G, leaves invariant. Choosing 

G = exp(e(t,P)), e(t,P) := J^tnP", (4.67) 

n>l 

with a constant N x N matrix P, and using ^(t, P)^[x] = C(t, P) + ln{lN - AP) with the iV x iV unit 
matrix T^, the heat hierarchy equations are mapped to 

{X - x_[^j)ix-' -p) = Xt, , (y - y^[x])ix-' -p)= yt, , (4.68) 

and (14.631 ) is converted into 

Xt^ +XP = RX + Qy . (4.69) 

Assuming that 

^(t,e)=A'(0)(T) + O(e), Kt, e) = 3^(o)(T) + 0(e) , (4.70) 

i? = i?(0) + O(e), Q = Q(0)+O(e), (4.71) 
and P independent of e, then we obtain from (14.691 ) 

e;f(o),Ti + ^(0)^ = ^{0)'^'{0) + Q(0)3^{0) + Oie) , (4.72) 

and from (14.681 ) 

e D{X)X^o) (A-^ -P) = eX(^o),T, + 0{e^) , (4.73) 
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together with the same equation for 3^(o) ■ After dividing the last equation by e, these equations have the 
dispersionless limits 

Rio)Xio)+Qio)ym=^io)P, (4.74) 

respectively 

Z)(A)A'(o) (A-i -P) = A'(o),T, , ^(A)3^(o) (A-^ -P)= 3^(o),Ti , (4.75) 

which is 

■;^^{o),T„+i = li^{o),T„P, ;^^^(o),T„+i = ^y(o),T„P, (4.76) 

(n = 1, 2, . . .), or in terms of the variables (13.18b . 

'^'(0),x„+, = X^oicc^P , 3^(o),,„^, = 3^(o),.„P , n = 1, 2, . . . . (4.77) 

Under the stated conditions, we have an expansion 

(A(t,e) = $(xi,X2,...) + 0(e), (4.78) 

which determines an exact solution $ of the dispersionless limit of the ncpKP hierarchy, i.e. the pdCM 
hierarchy ( 13.191 ). Proposition [T] in the following section confirms this directly, i.e. without reference to 
the scaling limit procedure applied to the ncpKP hierarchy and the above theorem. 

5 Exact solutions of the pdCM hierarchy 

The main result of the preceding section will be formulated in the next proposition, and we provide a 
direct proof. It will then be further elaborated and applied in order to construct some classes of exact 
solutions of the {su{m)) pdCM hierarchy. In this section, symbols like X and Q, for example, corre- 
spond to ^"(0) ^rid Q{0) in the preceding sections. Since now we resolve our considerations from the 
dispersionless limit procedure, there is no need to carry these indices with us any more. In fact, this 
section can be accessed almost completely without reference to the previous ones. 

Proposition 1 Let X be an invertible N x N and y an M x N matrix such that 

RX + Qy = XP (5.79) 

and 

= X^.P"" , = J^xiP" , n = 1, 2, . . . , (5.80) 

with constant matrices P, R of size N x N, and Q of size N x M. Then 

$ = yx~^ (5.81) 

solves the pdCM hierarchy 

^x„+UXrr.-^x^+^,x„ = [^x„,^xjQ 771,71 = 1,2,... (5.82) 

(which is ( 15. 79D with Q(o) replaced by Q). 
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Proof: (15.801 ) is equivalent to 

^x.„+, = x.„p , y.„^, = y.„p , n = 1, 2, . . . . 

In terms of the maps d, d, defined in section [331 this can be expressed as 

dA: = dXP, dy = dyp. 

Hence 

(d$) XP + ^dXP = dyP = dy = (d^>) X + ^dX = (d$) X + ^dXP, 

and thu^ 

d$ = (d$) W , 

where 

W :=XPX-^ =Q(^ + R^ 
using (I5.79I ). Since d and d satisfy (I3.44I ). and since Q and R are constant, we obtain 

dd$ = -dd$ = (d$) A dVF = (d$) A Q d$ , 

which is the hierarchy (I5.82I ). □ 

The next result shows how to obtain via proposition [T] solutions of the pdCM hierarchy in the algebra 
of m X m matrices with the usual matrix product (i.e. without the modification by a matrix Q different 
from the unit matrix). If ti{R) = tr(P), these solutions have values in sl{m, C). 

Proposition 2 Let U,V be N x m matrices and 

ip = U^<^V, Q = VU\ (5.83) 

where <!> = yX'^ with X, y solving (ETP]). Then 

tr{ip) = tr(P) - tr(i?) . (5.84) 

Under the conditions of propositionU] and if U and V are constant, ip solves the pdCM hierarchy dJ. 221 ). 

Proof: We have 

tr(99) = ti{u^yx-^v) = iT{vu^yx-^) = iT{QyX~^) . 

Using (15.791 ). this can be rewritten as 

ti{ip) = iv{XPX-^ - R) , 

which is (15.841 ). The last statement of the proposition is easily verified (see also section |3]). □ 
It is helpful to extend (15.791 ) to 

HZ = ZP, (5.85) 
^We note that this equation can be written as d$ — (d<I>)(5<I> — dO with O := "1>-R, which is l l3.16t . 
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where 

with constant matrices L, S. Indeed, the upper component of (15.851) reproduces (15.791) . But now we 
have an additional equation, namely SX + Ly = yP (which together with (15.791 ) implies the algebraic 
Riccati equation S + L<I> — — = for Although the latter appears to impose an unnec- 

essary restriction, it will be helpful in order to determine interesting classes of exact solutions. The two 
equations (15.801 ) can be combined into 

^x„+i =^xiP" n = 1,2,... . (5.87) 

Obviously, a transformation 

Z = TZ' , H = TH'T-^ , (5.88) 

with a constant matrix 

I''), (5.89) 

i 21 i 22 / 

preserves the form of the equations (15.851 ) and (15.871 ) with the same P. Consequently, if Z' solves (15.851 ) 
and (15.871 ) with H', and hence <!>' = y'X'~^ solves the pdCM hierarchy with Q', then Z solves the 
corresponding equations with H, and according to proposition [T] 

$ = yx-^ = {T21 + T22 $')(rii + ri2 «>T^ (5.90) 



solves the pdCM hierarchy with Qu Such a transformation thus relates solutions of different versions of 
the pdCM hierarchy, i.e. with different Q (which means different products). Since Q and Q' may have 
different rank, via (15.831 ) one obtains corresponding solutions of a pdCM hierarchy in a different matrix 
algebra. An extreme case is Q' = 0. Then the hierarchy (15.821 ) reduces to the system of linear equations 

X ^1 - X ^1 = m,n=l,2, .... (5.91) 

The above observation now suggests to first construct a solution of these linear equations, and then 
use such a transformation (as a "dressing transformation") to generate a solution of a nonlinear hierarchy. 

Proposition 3 Let P,L,R be constant N x N matrices. Let X' , y' solve ( 15.801 ) (which is ( 15.871) ) and 
Then 

<^ = y {X' - K yy^ (5.93) 

with any constant N x M matrix K, provided that the inverse in ( 15.951 ) exists, solves the pdCM hierarchy 
(ES21) with 

Q = RK - KL . (5.94) 



"'We note that the transformation ( 13.23b corresponds to the block-diagonal choice F — diag(gi2, <7i~^)- Such a transformation 
does not change a solution $ in an essential way. 

^If X' is invertible, then = y'X'~^ solves the linear hierarchy J5.91t . This follows from proposition [T] since Q' — 0. 
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Proof: Choosing in (15.881 ) the transformation matrix 
where In is the N x N unit matrix, we have 



R RK-KL 
L 



and hence Q = RK — KL. Since Z = TZ' again satisfies (15.851 ) and (15. 871 ). proposition [U tells us that 
$ given by ( |S!90l) . which is ( |S!93l ). solves (|5^ with Q given by ( |5^ . □ 

A special case of proposition [3] is formulated next. This will turn out to be particularly useful in the 
following. 

Corollary 1 Let P, K be constant N x N matrices, and X' an N x N matrix solution of 

K„+,=K,P'' n = l,2,..., (5.95) 

such that 

[P, X']=0. (5.96) 

Then 

^ = {X' - K)-^ , (5.97) 
provided that the inverse exists, solves the pdCM hierarchy with Q given by 

Q=[P,K]. (5.98) 
If moreover ^5.83^ holds, then ip solves the pdCM hierarchy l \3.22i in sl{m, C). 

Proof: We check that the assumptions of this coroUary constitute a special case of those of proposition [3] 
(15.851 ) decomposes into 

RX' = X'P , Ly' = y'p . 

Choosing 

R = L = p, y' = iN, 

this reduces to (|5^ . and (15^ reduces to (|5^ . Since ^'-^ = X' , ( |5^ becomes (|5^ . and (|5^ 
becomes (15.981 ). As a consequence of R = P and proposition [2l ip has vanishing trace, hence takes 
values in sl{m, C). □ 

Example. Let us choose 

with real constants pi ^ p2 and real functions /j. Then (15.961 ) holds and (15.951 ) requires that the function 
fi depends on the variables xi,X2, ■ ■ ■ only through the combination oji = Y^n>iP7~^^n- (15.981 ) is 
solved by 
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A diagonal part of K can be absorbed in (15.971) by redefinition of /i, /2. We obtain 

'^ = ^( r ^^^T^"), (5.101) 

V \ [Pl- P2) fl J 

where 

V = hf2 + {pi-P2)-\ (5.102) 
and then the following solution of the sl{2, M) pdCM hierarchy: 



1 / {P2-Pi) f: 



The corresponding conserved density £ is given by 

C_ 1+P1P2 rf/l df2 ^ ^^^^^ 

(/1/2 + {Pi -P2)~''^Y duJi duJ2 ' 
which can take both signs, depending on the values of the parameters. Choosing 

fi = exp{qi uJi) + Ci i = l,2, (5.105) 

with non-negative constants q and real constants q-i ^ 0, the solution is regular (for all xi, 2:2, • • •). For 
positive Cj, <S is exponentially localized, a sort of soli ton. The first derivatives of the components of ip are 
not localized, however. If c\ or C2 tends to zero, it stretches into a half-infinitely extended "line soliton", 
the location of which is determined by q\ uj\ + (72 1^2 = 0. □ 

As pointed out in section 13.11 the case where 99 given by (15.831 ) has values in the Lie algebra of a 
unitary group is distinguished by the fact that there is a non-negative "energy" functional, with density 
given by £ defined in (13.321 ). We will therefore concentrate on this case in the following. We further 
restrict our considerations to the case M = N , hence 3^, are all A'^ x matrices. Let U and V be 
constant N x m matrices. If <I> has the property 

$t = T^T'^ (5.106) 

with a constant invertible N x N matrix T which is anti-Hermitian, i.e. = —T, then by setting 

U = TV (5.107) 

we achieve that ip = U^^V is anti-Hermitian, i.e. 

ip^ = -ip . (5.108) 

As a consequence of these conditions, we have 

ip = -V^T^V, (5.109) 

and 

Q = VU^ = -VV^T , (5.110) 

which has the property 

Q^ = -TQT-^. (5.111) 
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We note that V ^ V a, with a constant unitary m x m matrix a, leaves Q invariant and induces a gauge 
transformation (p ip a. This can be used to reduce the freedom in the choice of V. 

In the following we address exact solutions of the su{m) pdCM hierarchy by using the recipe of 
corollary [T] Accordingly we should arrange that the solution X' of the linear hierarchy (15.951 ) satisfies 

X'"' = TX'T-'^ . (5.112) 

If also 

K^ = TKT-\ (5.113) 

then <I> given by (15.971 ) satisfies the same relation, i.e. (15.106b . As a further consequence, (15.1091) is then 
anti-Hermitian. 

Together with (15.1 121) . (15.961 ) implies [T^^ P^T, X'] = 0, which is identically satisfied as a conse- 
quence of ( 15.961 ) if P has the property 

P^ =TPT-^ . (5.114) 

We note that ( 15.981 ) is consistent with (15.1111 ). (15.1 131 ) and (15.1141) . Basically the problem of constructing 
solutions of (13.221 ) in su{m) (on the basis of corollary [T]) is reduced to the problem of satisfying the 
algebraic equation (15.981 ) with Q given by (15.1101 ). We summarize our results. 

Proposition 4 Let (P, X' , T, V) be data consisting of a constant N x N matrix P, an N x N matrix 
X', which solves ( 15.951 ) and ( 15.961 ). a constant anti-Hermitian N x N matrix T, and a constant N x m 
matrix V. Furthermore, let ( 15.7721 ) and l \5.114i be satisfied, Q be defined by ( 15.7701 ). and suppose that 
a solution K of ( 15. 981 ) and l \5.113\l exists. Then ip = —V^T^V, with ^ given by ( 15.971 ), is a solution of 
the pdCM hierarchy l \3.22i in the Lie algebra su{m). □ 

By application of proposition |4l some classes of exact solutions of the su{m) pdCM hierarchy will 
be derived in the following subsections. Examples are worked out for the su{2) case. Corresponding 
plots are restricted to the three variables entering the first hierarchy equation, and we will always use the 
coordinates t, x, y related to the variables xi, X2, X3 by the transformation (13.271 ). This is mainly done in 
order to ease a comparison with solutions of Ward's modified chiral model (cf. section [l!2l ). 



5.1 A class of solutions of the su{m) pdCM hierarchy 

Assuming that P is diagonal, i.e. 

P = diag{pi,...,pN), (5.115) 
with complex constants pi / pj for i / j, ( 15.961 ) requires X' to be diagonal. Writing 

;t" = diag(/i,...,/;v), (5.116) 
where the entries are functions of xi, 3:2, • • ., (15.951 ) becomes 

kx^=Pr^kx, j = l,...,N, n = l,2,.... (5.117) 
This is solved if fj is a holomorphi^ function of 

ivj ■.= ^XnP]-^ (5.118) 

n>l 

^More generally, the function fj is allowed to have singularities in the complex a;-, -plane, but we will not consider such 
solutions in this work. See also e.g. [44] in the case of Ward's chiral model. 
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(with the same j both for the function and its argument). In particular, fj depends on the variables 
xi,X2, ■ ■ ■ only through the combination (15.1 181) . The condition (15.1121 ) with an invertible anti-Hermitian 
matrix T imposes restrictions on the set of functions {fj}j=i,...,N, see section IF. 1.1 1 

According to corollary [T] ^ given by ( 15.971 ) solves the pdCM hierarchy with Q given by ( 15.981 ). 
which imphes Qa = 0,i = 1, . . . ,N, and 



Qi 



m N 

EE 



(5.119) 



where we took (15.1 101) into account. ( 15.971 ) shows that a diagonal part of K can be absorbed by redefi- 
nition of the functions fj. Hence it is no restriction to assume that Ka = for i = 1, . . . , A^. (15.1131 ) is 
then satisfied as a consequence of (15.1 141) . Now ( 15. 97b can be expressed as 



and ip = -V'^T^V solves the su{m) pdCM hierarchy (13.221) . 



(5.120) 



5.1.1 Some regular and localized solutions of the su{2) pdCM hierarchy 

Choosing even and for T the following block-diagonal form, 



/ -1 
1 



V 



-1 

1 / 



the conditions (15.1121) and (15.1141) read 



72(^2) = • • • , /Af(w7v) = fN-l{uJN-l)* : 

P2 =p*i,...,PN =P*N~1 ■ 

We also have uj2 = i^l, • • • , wat = 

Example 1. Let N = 2 and V = T2, which leads to Q = —T. Then we obtain 



P ( -i 

l + /32|/(a;)|2 V i 



(5.121) 



(5.122) 
(5.123) 



(5.124) 



where i = (3 = 2Q{p), p = pi, and / = /i is an arbitrary holomorphic function of a; = 

Yln>iP"'~^^n- This solution is regular for all 2:2, • • •■ The corresponding "energy density" is given 
by " 



|2\-2 



df 



djjj 



(5.125) 



Choosing for / a non-constant polynomial in cj, the solution is rational and localized, thus a field con- 
figuration that is often called a "lump". The shape of 8 depends on the degree of the polynomial and in 
particular on its zeros. 
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Figure 1: Plots of £ for stationary lump solutions of the su{2) pdCM equation, according to example 1 of 
section[5TI] Here we chose p = i and f{uj) =uj/2 (left), /(w) = (middle) and /(w) = (w^-l) /8 
(right). 



Let X4, X5, . . . = for the moment, so we concentrate on the first hierarchy equation. Applying the 
coordinate transformation (I3.27I ). we have 

uj = ^{t-x + 2py + p'^{t + x)) . (5.126) 

We note that this becomes t-independent if p = iti (i.e. (3 = ±2), in which case lo = —x ± i y, and the 
solution ip is stationary For f{uj) = + c, we obtain a simple lump. For example, choosing p = i and 

f{uj) = uj/2, we have 

» 

(5.127) 



(l+x2+y2)2 ' 

see figure [U c 7^ causes a displacement of the lump in the xy-plane. 

For f{uj) = q {uj — c)"', n > 1, with a zero of nth order, £ is bowl-shaped. In particular, ifp = i and 
f{uj) = tij^/2, we have 

_ 32 jx' + y^) 

^" (l + (x2+y2)2)2' ^^-^^^^ 

which is shown in figure [T] (second plot). The third plot in figure [T] displays another example. 

Configurations with M lumps are obtained by choosing / as a product of (powers of) factors lo — Ci, 
i = 1, . . . , M, with pairwise different complex constants Cj. □ 

Example 2. Let iV = 4 and 

y={ll)- (5.129) 

Then has the following components, 

^11 = -^22 = ^ [/3Mahlh2 - a*hh*2) - i + /?2)|6|^ + /3i|a/i2|' + /32|a/ii|']) , 
^12 = -<^2i = ^ /^i + a*/32h2) + a\h\^/32h2 + a*/3i/ii|/i2p) , (5.130) 

where 

Pi = 2 9=(pi) , a = pi-p2, b = pi-p2, 
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Figure 2: Plots of £ att = —2, 0, 2 for a three lump solution of the su{2) pdCM equation, according to 
example 2 of section U. 1.11 Here we chose pi = i,p2 = 2i, and /i(<jJi) = cof + 2, 72(1^2) = W2/4. Due 
to the special choice of pi, a pair of lumps is stationary. The positions of the latter are given by the zeros 
of fi{uJi) = — y^ + 2 — 2ixy, which are located at (x, y) = (0, ±\/2)- The position of the third lump 
corresponds to the zero of /2(w2)> which is given by (x, y) = (— 3t/5, 0). Choosing —2 instead of +2 in 
fi{uJi), all three lumps are located on the line y = 0, and the third lump moves through both members 
of the pair (which then reside at x = ±\/2)- 



hi 



■aPifi, h2 = a*(32f2, 

(|6|2 + |/li|2)(|6|2 + |/x2|2)+/3i/32|/ii-/i2 



(5.131) 



This solution is regular since V is positive (note that \b\ > since pi / p2, |6p > —Pi/32, and use 
+ |^2p > 1^1 — ^2p)- Figured shows an example. For generic parameter values, plots of £ show 
lumps with apparently trivial interaction. But if pi , p2 are close to the values ±1 (that correspond to the 
stationary single lump solutions), a non-trivial interaction is observed in a compact space region, see 
figure m The scalar KP-I equation possesses solutions with the same behaviour [85]. Moreover, also 
dipolar vortices (modons) of a barotropic equation [86] and BPS monopoles [87] show such a behaviour 
in head-on collisions. □ 



5.2 Another class of solutions of the su{m) pdCM hierarchy 

Let N be even. We introduce the commuting matrices 



fi hi 
// 



Pi 



PI 1 
PI 



(5.132) 



with pairwise different complex constants pj, and functions fi,hj, I = 1, . . . , N/2, and construct in 
terms of them the block-diagonal matrices 



V 



X2 



p 



/A 



V 



Pn/2 I 



which then obviously also commute. Now (15.951 ) becomes 

= vT^fl,x^ , hl,x„ = p'^~%,xr +{n- l)p"-V7,xi , 

where 1 = 1,..., N/2 and n = 1, 2, . . .. Writing 

df 



(5.133) 



hi = hi + 



dpi 



(5.134) 



(5.135) 
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Figure 3: Plots of £: at f = -110, -100, -99, -90, 0, 90, 99, 100, 110 for the solution in example 2 
of section [5. 1.11 with the data pi = —191/20, p2 = 21i/20, /i = 2\oJi, f2 = 2iuj2- Two lumps 
approach each other in x-direction, merge, move away from one another in y-direction up to some 
maximal distance, return to each other and merge again, and then separate in x-direction. 

the second equation is turned into = Pi~^hi,xi, by use of the first. Hence (15.951 ) is satisfied if, for 
/ = 1, . . . , N/2, fi and hi are holomorphic functions of 



UJI = '^XnP'^ ^ 



(5.136) 



n>l 



(which is (15.1 181 )). and in particular only depend on the variables xi,X2, ■ ■ ■ through this combination. 
In order to explore the consequences of ( 15.981 ). we write K and Q as N/2 x N/2 matrices, where the 
components Kjj, respectively Qij, are 2 x 2 matrices. 



Proposition 5 With the matrix P defined in ^5.133^ , and any Q, the solution of (i5.98i is given by 



for I ^ J, and 



where 



Kjj = _ [n2,Q/j] ^ [n2,[n2, Q/j]] 
PI - P.J {pi - pjY {pi - pjY 



[^2,Kjj]=Qjj J= l,...,A^/2, 



(5.137) 



(5.138) 



Ho 



1 




(5.139) 



20 



Dimakis and MuUer-Hoissen 



Proof: We write Pi =pil2 + Tl2. Then S5M\i . restricted 

(id H adn2)^/j = 

^ PI -PJ ' 

where oA^K = [U, K]. Now (15.1371) follows from 

(id + adn^) = V(-l)^(w-pj)-^ad^^, 

V pi-pj J 

since ad^j = 0. The diagonal components of ( I5.98I ) are adnj -f^jj = Qjj, which is (I5.138I ). □ 
Remark. In view of (15.971) . we may always assume that the two upper entries of Kjj vanish (since non- 
vanishing entries can be absorbed into X'). Using the matrix T given below in (I5.140I ). the condition 
(15.1 131 ) then implies k'!^j_-^ 2J-i ~ ''"2-f^2j,2j''"2> J = • • • > N/2, and this requires that Kjj can only 
have a non-zero entry in the lower left comer. As a consequence of (15. 1381 ). Qjj is then diagonal and 
has vanishing trace. □ 

A simple way of satisfying (15.1381 ) is to choose V such that the diagonal blocks Qjj vanish, and 
then set Kjj = 0, J = 1, . . . , N/2. This will be done in section [5".2. II 
It remains to satisfy the further anti-Hermiticity conditions. 



to components with / / J, takes the form 

^ Qij 
PI -pj' 



5.2.1 su{2) lumps with "anomalous" scattering 

Let N now be a multiple of 4. In analogy with (15.1211 ) we set 



/ 

T2 



\ 



-T2 





Then (15.1141 ) means pI = T2P2T2, 





P^ 

' ^ N/2-1 



-T2 

/ 



where T2 



1 

1 



(5.140) 



P2J = P2.J-I 



T2PN/2T2, which by use of (15.1321 ) amounts to 

J = l,...,iV/2 . (5.141) 



Since we address the case m = 2,V has to be chosen as an x 2 matrix, which we subdivide into 
2x2 blocks V/, / = 1, . . . , N/2. It follows that 



Q 



J.I 



V ^ if J IS 



r2 



odd 
even 



Thus, in order to achieve that Qjj = 0, we must arrange that 



V2J-1VI = 
Then Q has the following structure 



Q 





-V2V^ 







J = l,...,N/2 



-V2VI V2VI 



(5.142) 



(5.143) 



T2 



(5.144) 
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Example. The simplest case is = 4. Excluding degenerate cases, the two blocks Vi, V2 of V should 
both have rank 1. Hence Vi = viu\, V2 = V2u\ with vectors uj, vj, J = 1,2, satisfying u\u2 = 0. With 
a unitary transformation a we can achieve that the lower component of ui vanishes. It follows that the 
upper component of U2 also vanishes. By a redefinition of ^1,^2, we obtain u\ = (1, 0) and = (0, 1), 
and thus 

^ ( \] ( V11VI2 blip 

V 1^121 ^^12^^11 

Q2I = -V2{T2V2y = - I f2 ' , (5.145) 

V I ^22 1 ^^221^21 / 

with an obvious notation for the components of vi and f 2. We should exclude the case when the expres- 
sion (15.1371 ) for K reduces to the first term on the right hand side, since this leads back to the solution 
of example 1 in section [5TT] This case is ruled out if V12 or V22 is different from zero, which suggests to 
choose = (0, 1), J = 1, 2, and thus 



Then proposition |5] yields 



where 



-{Pi-Pj) ^ -2{pi-pj] ^ 
{Pi-Pj)~^ [pi-Pj] 



Kij = [ r :\u ir' J(U I (5.148) 



for I ^ J. This in turn allows to compute $ and then also 99. The anti-Hermiticity conditions are then 
satisfied by setting 

P2=PI, /2(W2) = /1(W1)*, h2{uJ2) = hi{uj{)\ (5.149) 

and we have a; := cji = uJ2- The result is 

ifii = -9^22 = (2 + + /3V + ' 

¥'i2 = -y^2i = ^(4i/-/3^-/3W) , (5.150) 

with (3 = 29(j)i), / = fi{io), h = hi given by (15.1351 ) in terms of / and h = hi{u}), and 

V = r'(i + aV?? + r'|2/ + i/3/ip . (5.151) 

The solution Lp is thus regular for any choice of pi, with non- vanishing imaginary part, and the holomor- 
phic functions /, h. An example with /i = is shown in figure |4l More interesting structures appear 
for non-constant h. Indeed, figure |5] shows two lumps that scatter at an angle of 90°. Choosing / linear 
in Lo and h proportional to w", we observe a vr/n scattering. Figures |6] and |7] show examples of 60°, 
respectively 45° scattering. □ 

Solutions with vr/n scattering have also been found in Ward's chiral model numerically [35, 37], 
and analytically as certain limits of families of non-interacting lumps [32,39,43,50]. Moreover, also 
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Figure 4: Plots of £^ at t = —50, 0, 50 for the solution of section [5^.2. ll with the data pi = i (i.e. /3 = 2), 

f = 2uja.ndh = 0. 




Figure 5: Plots of f at t = —5,0,5 for a 2-lump solution, exhibiting "scattering at right angle", see 
section IT.2. II Here we chose pi = i and f = —iLu/32, h = — a;^/4. 



the scalar KP equation (with positive dispersion, i.e. KP I) possesses solutions with this behaviour 
[88-93]. In fact, vr/n scattering in head-on collisions of soli ton-like objects is a familiar feature of 
many models (see [94-96], in particular). It occurs in dipolar vortex collisions [86,97-99], in 0(3) 
and CP^ models [37, 100-106], in Skyrme models [95, 106-109], for vortices of the AbeUan Higgs (or 
Ginzburg-Landau) model [106, 110-119], and BPS monopoles of a SU{m) Yang-Mills-Higgs system 
[87, 106, 120-122]. Another integrable system that possesses solutions with this behaviour is the Davey- 
Stewartson II equation [123, 124] (which can actually be obtained by a multiscale expansion from the 
KP equation [24, 125]). 

The fact that lumps can interact either trivially or non-trivially (in Ward's chiral model) has been 
attributed to the status of the internal degrees of freedom in the solutions [32]. But such an explanation 
appears not to be applicable to the case of the scalar KP equation. This requires further clarification. 



5.3 A further generalization 



In case of the solutions obtained in section [5^21 the matrix P consists of complex conjugate pairs of 2 x 2 
blocks of Jordan normal form. Of course, this can be generalized to Nj x Nj Jordan blocks 



f PI I • • • \ 

PI 1 ■•• i 



Pi 



\ 



PI 1 
PI j 



(5.152) 



Dispersionless limit of the ncpKP hierarchy 



23 




Figure 6: Plots of <S at t = —500,0,500 for a 3-lump configuration exhibiting 60° scattering. This is 
obtained from the solution of the example in section l5".2.1l with the data pi = i and f = —ilo, h = lo'^/S. 




Figure 7: Plots of f at f = —20, 0, 20 for a 4-lump configuration exhibiting 45° scattering. This is 
obtained from the solution of the example in section |5".2. II with the data pi = i and f = —ilo, h = 



and P can be chosen as a block-diagonal matrix with pairs of conjugate blocks of this form. For each 
pair (P/ , Pj ) in P, the matrix T should then have a corresponding block 



Ti 



V 1 



-1 \ 



(5.153) 



of size 2Nj x 2Nj, in order to achieve that (15.114b holds. 
Example. Let 



P 



(P 1 
p 1 

p 



p* 1 

p* 1 
Op*/ 



\ ( 



-1 \ 



V 1 



V 



( 



1 




V 



\ 







1 / 



(5.154) 
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Figure 8: Plots of <S at t = —20, 0, 20 for the solution of the example in section [531 with the data p = i, 
f = —luj, h = w^/S and 5 = 0. At t = we have cut off an extremely large lump in the center. 



Then X' must have the form 



X' 



f f h g 
/ h 
f 



f* h* g* 
f* h* 

r J 



(5.155) 



and (15.801 ) becomes 



9xn 



X, + (n - l)p" 



JXl 5 



Writing 



h = h 



dl 
dp 



9 = 9 + 



dh 1 



(5.156) 



(5.157) 



dp 2 dp'^ ' 

with functions f,g,h, it follows that these equations are satisfied if the latter are arbitrary holomorphic 
functions of a; = Yln>i P^^n- Furthermore, we find that 



K 





-"■12 



K12 





with K12 





(5.158) 



and /? 



: 29 (p), solves {P,K\ 
det(A"-i^) = /? 



= Q with Q = — VV^T. The resulting class of solutions is regular since 

+p-'\2f + (/ + iphf + /3'f~g\' + /3-i2|6/ + p^4ih- P~g)\' . (5.159) 



Now we have three arbitrary holomorphic functions at our disposal, so this class exhibits quite a variety 
of different structures. Figures [8] and |9] show some examples. If h = g = 0, the typical behaviour is 
similar to the one shown in figure m □ 

Comparing the data that determine the class of solutions in the example in section 15.2.11 based on a 
conjugate pair of 2 x 2 Jordan normal form matrices Pj, with those of the last example, which is based 
on a conjugate pair of 3 x 3 Jordan normal form matrices, there is an obvious generalization to the case 
of conjugate pairs of larger Jordan normal form matrices Pj. We note in particular that proposition [5]can 
be generalized. Since the solutions turned out to be automatically regular in the 2x2 and 3x3 case, it 
may well be that this holds in general. But a proof of this conjecture is out of reach so far. 



Dispersionless limit of the ncpKP hierarchy 



25 




Figure 9: Plots of <S at t = —1000, 0, 1000 for the solution of the example in section [531 with the data 
p = i, f = —iuj/AO, h = u!^/W and (7 = 0. At t = we have cut off the lumps beyond a certain height. 



5.4 Superposing solutions 

The data that determine solutions of the su{m) pdCM hierarchy on the basis of proposition |4] are given 
by a set of matrices (P, X' , T, V). Let two such sets be given, (Pj, A'/, Tj, Vi), i = 1,2, with associated 
matrices Qi = -ViV^Ti, Ki (as solutions of (15.981) ). and <I>j given by (15.971) . Pi,Xi,Ti are Ni x Ni 
matrices and 14 is an A'^j x m matrix. We can combine them into the larger matrices 

X{ 



P 

T : 



Pi 

P2 

Ti 

T2 



X' 

V -- 



X!, 
Vi 
V2 



(5.160) 



Obviously, (P, X' , T, V) again satisfies ^5^ . dlH, ^5TY2\ . (l5Tl4l) and Tt = -T. KTM becomes 



Q 



Qi Q12 
Q21 Q2 



and iSWi with 



yields the equations 



K 



Q12 = -VxVi T2 , Q21 = -V2V; Ti , 



Ki K12 
K21 K2 



(5.161) 



(5.162) 



P1K12 - K12P2 = Q 



12 



P2K21 - K21P1 = Q 



21 



(5.163) 
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The off-diagonal blocks of K are a source of complexity and non-triviaUty of the resulting superposition. 
(15.1131 ) then determines K21 in terms of K12 (or vice versa), 

i^2i = T2-i/4ri . (5.164) 

As a consequence, the second of equations (15.1631) follows from the first. If we find a solutioij^ K12 of 
the remaining equation, then we obtain 

" = [ <^.K2,^. ^2 )[ A,- ) ' ^'-'^'^ 

where 

Ai = In, - Ki2^2K2i^i , A2 = 1n2 - K2i^iKi2^2 , (5.166) 

and Lp given by ( 15.1091 ) solves the su{m) pdCM hierarchy, provided that the inverses of Ai and A2 exist. 
If the two matrices <I>j are regular (and thus also the corresponding solutions Lpi), then <I> and thus also 
is regular if and only if det(Ai) det(A2) 7^ (for all values of xi, X2, • • •). Since det(Ai) = det(j42) 
by an application of Sylvester's determinant theorem, this reduces to the condition 

det(Ai) / . (5.167) 

We note also that det(Ai) is real since 



det(Ai)* = det(4) = det(TiVi - kI-^<^\k\^^\) 



det(jjvi - TiKi2^2K2i^iT{^) = det(^i) . (5.168) 



Example. We choose 



Ti 



( Pi 


V 



V 1 



1 

Pi 



Pi 





1 

p\ J 



( h hi + dfi/dpi 
/i 



V 



f! K + idfi/dpiY 
A* 



/o o\ 
1 

V 1 y 



where /i, hi are arbitrary holomorphic functions of uji (with lui defined in (15.1 181 )). and 



P2 



P2 







P2 



f2 







/2 



Vo 



1 
1 



(5.169) 



(5.170) 



with an arbitrary holomorphic function /2 of 102- Thus we superpose data corresponding to a regular 
solution of the kind treated in the example in section 15.2.11 and data corresponding to a regular solution 
as given in example 1 of section |5?T] In the following we assume that 



Pl7^P2, P\+P2- 



(5.171) 



'choosing Vi and V2 such that V\V^ = 0, we have Q\2 — and l l5.163t is solved by 7^12 = 0. It follows that ip is simply 



the sum of the solutions cpi and ip2- But V1V2 
be degenerate, i.e. cannot have full rank. 



also implies that 951^2 = 0, hence both constituent solutions ipi, ip2 must 
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Together with the conditions pi ^ p*, i = 1,2, which the data of the components have to satisfy, this 
means that the constants pi and their complex conjugates are pairwise different. The second condition 
in (15.1711 ) is in fact needed for the matrix K to exist. K has the form (15.1621 ). where Ki is given by the 
4x4 matrix K in (15.1471 ) with the pair {pi,Pi). Furthermore, 



K2 





1 



P2-P2 





21 







1 



1 



Pi -Pi 



(pi-p|)2 



Pj-Pa ipl-P2r 




(5.172) 



and Ki2 is then determined by (15.1641 ). With some efforts the expression for det(j4i) can be brought 
into the form 



det(^i) 



a\-\l + f3l\f2\')-' (kP + (1 + ^flhm-' [{\a\'[\w\' + (1 + PtlM'f ] 
-P,I32 [ |ap(l + /3f + + \aw + ib*l3fh\'' ])' 
{a*)\\w\^ + {l + Pt\hm(32f2 



+ 



+Pi^2{b'w + ia*f(3fffw* - 2ia*b0tfia+Pt\fin 



where 



(3i = 2Q{pi), a = pi-p2, b = pi-p2, w = (3fhi - 2i(3ffi 



(5.173) 



(5.174) 



The regularity condition (15.1671 ) turns out to be automatically satisfied. This is seen as follows. First we 
note that 



\a\' = \b\' + P1P2 (|a|2 + > (5,P2 {\af + \bf) 
as a consequence of the first of the inequalities (15.1711) . and thus 



ar (1 + Ptifii'f > Pi/32 iH' + mil + ptih 



(5.175) 



(5.176) 



Using |ap > /3i/?2, this leads to 



w 



+ + > ^,p,{\af{l + pf\f,\^) + \b\' + \aw\' + \ib*Plh\^) 

> Pi(32 {\af{l+pf\fi\^) + \b\^ + \aw + ib*(3!fif) , (5.177) 

which implies det(Ai) > 0. 

Figure [To] shows plots of £ at consecutive times, for a special choice of the data. □ 

In the last example the regularity of the superposition turned out to be a consequence of the "regular 
data" we started with. But this example also demonstrates that it is quite difficult in general to evaluate 
the regularity condition (15.1671 ). We note that also the cases treated in sections [5 . 1 . 1 1 and [5.2.1 I may be re- 
garded as special cases of "superpositions" as formulated above. In particular, example 2 of section [5". 1.1 1 
provides us with another example where the superposition of regular data turned out to be regular again. 
It is unlikely that this is a special feature of our particular examples. But in order to tackle a general 
proof, we probably need different methods. 



6 Conclusions 



We summarize the relations between integrable systems and their hierarchies considered in this work in 
the following diagram. 
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Figure 10: Plots of iS at t = —2, 0, 2 for a superposition of a 2-lump configuration, with "anomalous 
scattering", and a single lump (which is at the top of the left plot and at the bottom of the right plot), 
according to the example of section l54l Here we chose pi = i, p2 = — 3i/8, /i = — ia;i/32, /2 = 24^2, 
and hi = —uj"^ /A. 



dispersionless limit 
mKPg Wardg 

Miura transf . | | pseudo-duality 

t dispersionless limit t 
pKPq pdCMQ 

rank((5) = 1 | reahty cond. rank(Q) = | su{m) cond. 

scalar pKP su{m) pdCM 

Here pKPg and pdCMg stand, respectively, for the pKP and pdCM hierarchy in the matrix algebra 
GL{M X N, C) with product modified by a constant matrix Q (see (14.621) ). pdCMg is related by pseudo- 
duality (see (13.361) ) to the hierarchy Wardg of Ward's model with dependent variable in GL{M x N, C) 
(and product modified by Q). If rank((5) = 1, solutions of pKP^ are mapped to solutions of the 
scalar pKP hierarchy, an additional condition ensures that the resulting solution is real. Analogously, if 
rank((5) = m and an su{m) condition holds, solutions of the pdCMg hierarchy are mapped to solutions 
of the su{m) pdCM hierarchy (which is pseudo-dual to the hierarchy associated with Ward's modified 
SU{m) chiral model). Concerning the Miura transformation between the (matrix) pKP hierarchy and 
the modified KP (mKP) hierarchy, and its dispersionless limit (see the dashed arrows in the diagram), 
see [76]. The relations provided by the dispersionless scaling limits in the diagram have actually been 
anticipated in [76] (see the remark in section 4 therein). 

In the present work we demonstrated how the dispersionless scaling can be used to transfer a method 
of constructing exact solutions from the (matrix or "noncommutative") pKP hierarchy to the pdCM 
hierarchy. Indeed, proposition [T] is an analog of theorem 4.1 in [75] (which we recalled as theorem [T]). 
We showed that large classes of exact solutions of the pdCM hierarchy can be obtained with its help. 
In particular, we presented examples of various multiple lump configurations of the su{m) pdCM. The 
general result formulated in proposition [T] is a source of even more classes of exact solutions. 

Our method to generate exact solutions of the su{m) pdCM hierarchy is based on quite simple for- 
mulae and quickly produces interesting solutions (like lumps with "anomalous scattering"). But a more 
systematic treatment, in particular of multi-lump solutions, requires deeper methods (of matrix calcu- 
lus), and further insights are needed as to how the a priori given plethora of parameters can efficiently be 
reduced. It would also be of interest to compare this method with an inverse scattering approach. 

Solutions of the su{m) pdCM hierarchy can also be obtained from solutions of Ward's chiral model 
hierarchy by integrating (13.361 ) (or equivalently (13.511 )). In any case, one should expect an analogous 
structure of localized solutions, and this indeed turns out to be the case in examples. A deviation in 
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the corresponding plots is caused by the fact that the "energy" expression for the pdCM differs from 
the energy of the Ward model by a term that causes an asymmetry in the x-direction, see the remark in 
section 13.21 and also the appendixlf] Our method to generate solutions of the pdCM hierarchy seems to 
be quite different from the methods that were used to construct solutions of Ward's model. In particular, 
in the latter model solutions with "anomalous scattering" have been obtained by taking suitable limits 
of families of non-interacting lump solutions. In our approach, corresponding solutions of the pdCM 
hierarchy are directly given by matrix data involving Jordan blocks. Moreover, we have seen that even 
the simple multiple lump solutions of section 15.1.11 can exhibit an anomalous behaviour within some 
compact space region (see figureO, whereas asymptotically (i.e. compared at large enough negative and 
positive times) no deflection is observed. We noted that this has a KP-I counterpart [85] and also analogs 
in some other systems [86, 87]. 

The fact that the dispersionless scaling limit of matrix pKP (respectively mKP) is a simple reduction 
of a potential version of the (4-dimensional) self-dual Yang-Mills equation raises the question whether 
there is a (4-dimensional) integrable system that has the full self-dual Yang-Mills equation as a disper- 
sionless limit and that admits a reduction to matrix KP (respectively mKP). 

Appendix: The "anti-pdCM hierarchy" 

Writing (l330l) as Jd{J-^dJ)J-^ = 0, using the Leibniz rule and (13.441 ). we obtain the equivalent form 

d((dJ)J-i) = (1.178) 

of the hierarchy associated with Ward's chiral model, where d and d exchanged their roles. This is 
integrated by introducing a potential ip such that 

{dJ)J-^=dif. (1.179) 

Rewriting the last equation in the form dJ = (d(p)J, we obtain the integrability condition 

dd(^ = A . (1.180) 

In components, this becomes 

'fx„X,r.+ l - 'fx^X,, + l = ['fx„ + -,,'fx,r,+ l] 171, n= 1,2,... , (1.181) 

which we refer to as the anti-pdCM hierarchy. For m = 1 and n = 2 we have 

<i^xia;3 — 'f>X2X2 = ['^X2^'^Xi] ■ (1.182) 

In terms of the coordinates given by (I3.27I ). it takes the form 

'^tt - ^xx - ^yy + {^t + (^x, (^y] = . (1.183) 

Since this equation is obtained from (13.281) by x —x, so are its Lagrangian £ and energy-momentum 
tensor T^i, from those in section ITT] For (p in su{m), 

£ = f\ + f\ = -hi^i^t + ^x? + (1.184) 

is then a non-negative conserved density. 

'^We note that x and y have to be exchanged for comparing our formulae with those in the literature on the Ward model. 
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Associated with any solution J of Ward's chiral model hierarchy via (13.511 ) and (I1.179I ). there are 
solutions if and ip of the pdCM hierarchy and the anti-pdCM hierarchy, respectively. Using JtJ^^ — 
JxJ^^ = ify and JyJ~^ = (pt + (px, which follow from (11.1791) . we find that 

£ = <?Ward + tr(j-i Jt J-^Jx) . (1.185) 
Combining this with (13.411) . leads to 

£Y^.rA = \{£ + £) . (1.186) 

The next result is an analogue of corollary [T] and leads to a class of solutions of the anti-pdCM 
hierarchy. 

Proposition Al Let (P, X') be data that determine via corollary\l\a solution $ (given by ( 15.971 )) of 
the pdCM hierarchy with Q = [P, K]. If P is invertible (as in all our examples in section\5^, then $ also 
solves 

dd$ = d$AQd^> (1.187) 

with Q = p-^QP-\ 

Proof: As a consequence of ( 15.951) . X := X' — K solves 

dX = dxp-^, 

and consequently ^ = X^^ satisfies 

d$ = (d$)VK-\ w = xpx-^. 
Now we note that K96\i and (15^981) imply [P~'^,X] = [K,p-^= p-'^QP-^. Hence 
^-1 ^ + QP~^) X~'^ = P-^ + P-^QP-^^ , 

and we obtain 

dd^> = -dd$ = (d$) A p-^Qp-^d<^ . 

□ 

If moreover the assumptions of proposition |4] are satisfied, then 

P-^QP-^ = -p-^VV^TP-^ = -P-^V{P-^V)^T , (1.188) 
and (11.1871 ) imphes that 

= -(P"V)^r$P"V (1.189) 

solves the anti-Hermitian anti-pdCM hierarchy (11.1801 ). The data {P, X' ,T,V) therefore determine a 
solution (15.1091 ) of the anti-Hermitian pdCM hierarchy and also a solution (11.1891 ) of the anti-Hermitian 
anti-pdCM hierarchy. 

Although elaboration of examples suggests that the pair {(p, (p) determined by the data (P, X' , T, V) 
indeed belongs to the same solution J of Ward's chiral model hierarchy (via (13.511 ) and ( 11.1791 )). we were 
not able so far to prove this. 
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